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EFFECTS OF BROWSING BY NATIVE UNGULATES ON THE SHRUBS IN BIG 
SAGEBRUSH COMMUNITIES IN YELLOWSTONE NATIONAL PARK 


Francis J. Singer! and Roy A. Renkin? 


ABSTRACT.—The effects of elk (Certus elaphus), pronghorn (Antilocapra americana), and mule deer (Odocoileus 
hemionus) browsing on shrubs in big sagebrush (Artemisia tridentata) communities were monitored over a 31-year peri- 
od in Yellowstone National Park. Ungulates were restricting Wyoming big sagebrush (spp. wyomingensis) heights, size, 
and recruitment on the lower-elevation stratum only, while no such suppression was observed on the high-elevation 
stratum. Parallel increases in mountain big sagebrush (spp. raseyana) densities and cover occurred over the study period 
on both browsed and unbrowsed sites at the higher-elevation stratum, although big sagebrush, green rabbitbrush 
(Chrysothamnus viscidiflorus), and horsebrush (Tetradymia canescens) were slightly taller and crown sizes were slightly 
larger on unbrowsed than browsed sites. Wyoming big sagebrush utilization (percent leader use) was eight times higher 
(x = 87 + 7.2% by pronghorns, mule deer, and elk) on the low-elevation winter range stratum (the Boundary Line Area 
[BLA] portion of the winter range), while mostly mountain big sagebrush with leader use averaged only 11 + 4.1% (near- 
ly all by elk) on the high-elevation range stratum. In addition, annual aboveground biomass production of big sagebrush 
did not differ between browsed and unbrowsed study sites on the high-elevation stratum of the winter range. Population 
turnover was higher on browsed big sagebrush at the high-elevation plots; seedling germination and survival rates were 
higher on browsed plots versus unbrowsed plots. No difference was observed in percent dieback of big sagebrush adult 
plants between browsed and unbrowsed plots at the higher stratum. Browsing did not influence the number of leaves or 
seedstalks per plant (P > .05), but leaves averaged 45% longer and seedstalks 42% longer on browsed big sagebrush. 
Ungulate browsing, however, apparently suppressed production, germination, and survival of Wyoming big sagebrush 
on the low-elevation stratum. Numbers of Wyoming big sagebrush declined 43% and cover declined 29%, 1957-1990, 
on browsed sites on the BLA. Annual biomass production on browsed sites at the low-elevation stratum was only 6-35% 
that of unbrowsed sites, and big sagebrush recruitment was Jess on browsed sites. Percent leader use of big sagebrush 
did not differ between the period of ungulate reductions, 1962-1969, and the 1980s on the Jower stratum (¥ = 87% 
leader use), but utilization was less on higher portions of the winter range during the period of elk reductions (¥ = 2%) 
than during the 1980s following cessation of elk controls (¥ = 11%). 
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Native populations of elk (Cervus elaphus), 
bison (Bison bison), and pronghorn (Antilocapra 
americana) were artificially reduced in Yellow- 
stone National Park (YNP), particularly from 
1942 through 1967 (Meagher 1973, Houston 
1982). Reductions were terminated in 1967 
when an experimental management program 
of natural regulation was initiated (Cole 1971, 
Houston 1976, 1982). Elk and other ungulate 
numbers tripled after cessation of controls, and 
concerns were expressed over high ungulate 
densities (Chase 1986, Kay 1991). Appropriate 
numbers of ungulates for the park are unknown 
since no similar control area exists where wolves 
(Canis lupus) are present and where ungulate 
migrations are completely unrestricted by 
humans (Cayot et al. 1979, Peek 1980). Ungulate 
densities are likely slightly above natural con- 
ditions, in that three independent computer 


models suggest 8-15% fewer elk and 10-25% 
fewer bison would occupy the system if wolves 
were recovered (Garton et al. 1990, Boyce 1993, 
Mack and Singer 1993). 

Early workers expressed concern about 
apparent overbrowsing and declines in big sage- 
brush (Artemisia tridentata) due to possible 
overabundant populations of elk and prong- 
horn. As early as the 1930s, Rush (1932) and 
Cahalane (1943) reported losses of big sage- 
brush over lower-elevation areas of the north- 
erm winter range. Rush (1932) reported that 
less-palatable rabbitbrushes (Chrysothamnus 
spp.) were increasing. Kittams (1950) concluded 
that big sagebrush numbers were declining at 
both lower and higher elevations of the north- 
erm winter range. He felt that physical distur- 
bances of big sagebrush by elk during cold 
periods (shattering and trampling) and an 
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absence of big sagebrush reproduction con- 
tributed to the decline. Declines in big sage- 
brush at the lower-elevation Boundary Line 
Area (BLA) were attributed by Kittams (1950) 
to excessive levels of browsing by pronghorn. 
Park management established a goal to reduce 
the pronghorn herd by 50% (Kittams 1959); by 
1969 pronghorn numbers were artificially re- 
duced from 600-S00 to less than 200 through 
a combination of artificial reductions and 
severe winters (Barmore 1980). 

Houston (1982) provided alternative inter- 
pretations concerning big sagebrush. He re- 
ported increases in big sagebrush numbers over 
all the northern winter range except the BLA 
near Gardiner, MT. where numbers declined. 
Houston (1982) compared photos taken during 
the 1860s to photos retaken in the 1970s. He 
attributed the increase in big sagebrush at 
higher elevations to fire suppression and the 
decline in big sagebrush in the BLA to a 
return to more natural conditions following 
the removal of intense grazing by livestock in 
the early 1930s when the area was added to 
the park. In 1986 the U.S. Congress directed 
the National Park Service (NPS) to conduct a 
study to evaluate whether native ungulates 
were overgrazing the northern winter range 
(Congressional Record 1986). 

Our objectives were to document trends in 
big sagebrush abundance on a series of per- 
manently marked plots from 1958 to 1990. 
Height, canopy size, twig lengths, and annual 
production of shrubs were compared between 
browsed and unbrowsed sites. 


STUDY AREA 


Shrub sampling was conducted on un- 
browsed (exclosed) and paired browsed sites at 
eight ungulate exclosures erected in 1958 and 
1962 on Yellowstone’s northern winter range. 
The eight exclosures, 2 ha in size, were locat- 
ed on gently rolling upland steppe ridge and 
the intervening swale habitats (Fig. 1). We 
divided the study area into a low-clevation 
stratum (the BLA of Houston 1982) with two 
exclosures, and a much larger, high-elevation 
stratum (n = 6 exclosures), based on large dif- 
ferences in ungulate species, elevation, snow- 
pack, precipitation, and big sagebrush sub- 
species (Fig. 1). 

Underlying soils are typic ealciborolls, aridic 
haploborolls, and aridie calciborolls (Lane 


GREAT BASIN NATURALIST 


[Volume 55 


1990). Precipitation averages 30 cm at the low- 
clevation exclosures and 55 cm at the high- 
stratum exclosures (Ifouston 1982, Despain 
1991). The northern winter range is approxi- 
mately 1100 km? and is located along the 
upper Yellowstone River drainage. Elevations 
range from 1500 m at the fow-stratum exclo- 
sures to about 2200 m at the highest exclosure. 
The northern winter range is lower, warmer, 
and drier than the remaining higher plateaus 
of YNP (Houston 1982). As a result, 80% of the 
ungulates in the park during winter are found 
on the northern winter range (Singer 1991). 

Dominant shrubs at all the study sites 
include two subspecies of big sagebrush at the 
high-elevation stratum, nearly all mountain big 
sagebrush (A. ¢. vaseyana) with some basin big 
sagebrush (Artemisia tridentata tridentata), while 
Wyoming big sagebrush (A. t. wyomingensis) 
occurs in the low-stratum study sites. The big 
sagebrush subspecies vary markedly in their 
site requirements, growth, and preferences by 
ungulates (Beetle 1960, Welch et al. 1981, Beetle 
and Johnson 1982, McArthur and Welch 
1982). Rubber rabbitbrush (Chrysothamnus 
nauseosus), green rabbitbrush (Chrysothamnus 
viseidiflorus), and horsebrush (Tetradymia 
canescens) are found at all study sites. Grayia 
spinosa and Atriplex eaneseens occur at the 
low-stratum study sites (Houston 1982), 
Dominant grasses are bluebunch wheatgrass 
(Pseudoroegneria spicata), tdaho fescue (Festuca 
idahoensis), janegrass (Koeleria pyramidata), 
bluegrasses (Poa compressa, PB. sandbergii, or P. 
pratensis), and thick-spike wheatgrass (A. 
dasystachyum; Houston 1982, Wambolt et al. 
1987, Despain 1991). 

Pronghorn and mule deer occupy only the 
low-elevation stratum of approximately 52 km? 
within Yellowstone National Park (Barmore 
1980, Houston 1982, Singer 1991). About one- 
half of the pronghorn population also sum- 
mers in the low-elevation stratum. The lower- 
elevation (about 1500 m) exclosures are locat- 
ed in typical, nearly snow-free, rolling xeric 
shrub and mixed grassland habitats. Elk occu- 
py both strata and the entire winter range of 
§10-1000 km2 (Houston 1982, Singer 1991). 
Only elk and bison winter near the high-eleva- 
tion stratum exclosure sites. Winter snow depths 
near the higher exclosures (1639-2200 m) are 
typically 0.4-0.6 m, which are excessive for 
pronghorn and deer. The abrupt elevation rise 
for Mt. Everts separates the high and low strata. 
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Fig. 1. Map of the northern Yellowstone ungnlate winter range and the high-elevation stratum (six exclosures) and 
low-elevation stratum (two exclosures) big sagebrush study sites. Pronghorn, mule deer, and the Wvoming subspecies of 


big sagebrush were found only at the low-elevation stratum. 


The bison winter range expanded from about 
130 km? in the 1960s in the higher stratum to 
about 460 km? in the late 1980s during a peri- 
od of bison population and range expansion 
(Meagher 1989, Singer and Norland 1995). 
Periodic bison use of the low-elevation stratum 
occurred following the population expansion 
in the late 1980s (Meagher 1989), but vegeta- 
tion measures reported here are nearly all prior 
to any bison use of the low-elevation stratum. 
Elk numbered about 8000 in 1958-1962 
when initial monitoring of the sagebrush belt 
transects began. Elk were subsequently further 
reduced by artificial controls to less than 5000 
in 1967 (Houston 1982). After cessation of con- 
trols, elk steadily increased, with counts rang- 
ing from 16,000 to 19,000 from 1982 to 1989 
(Singer et al. 1989, Singer 1991). Bison were also 
artificially controlled until 1967. After cessa- 
tion of controls, bison on the northern range 
increased from less than 100 in 1967 to 850 by 


1988 (Houston 1982, Meagher 1989). Prong- 
horn were artificially reduced from 600-800 to 
<200 (Barmore 1980), and pronghorn num- 
bers remained <200 until about 198]. During 
the 1980s—apparently due to milder win- 
ters—pronghorn increased to about 600 (Singer 
1991). Mule deer counts increased from 1000 
in 1985 to 2300 in 1988 over the entire deer 
winter range, the majority of which lies north of 
the park boundary (Singer 1991). Conversely, 
mule deer counts just within the park bound- 
aries declined from 230 in the 1960s (Barmore 
1980) to about 100 in 1988, in spite of the 
overall herd increase. 





METHODS 
Ungulate Densities and Diets 


Average ungulate densities near the exclo- 
sures were based on actual aerial counts made 
from fixed-wing aircraft as described in 
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Barmore (1980), Houston (1982), Meagher 
(1989), Singer (1991), and Singer and Norland 
(1995). Densities are uncorrected for visibility 
bias (Samuel ct al. 1985) and therefore repre- 
sent minimum average densities for the study 
periods—undoubtedly some animals were 
missed on the counts (Singer et al. 1989). Diets 
of all four ungulates found near the study sites 
on the northern winter range were estimated 
for each of three winters, December—March 
1985-1988, from microhistological analysis of 
fecal samples (Washington State University, 
Wildlife Habitat Laboratory, Pullman). Each 
sample was a composite of 5 g of fresh dung 
material from 6-12 dung piles. Aggregate 
average percentages are reported for signifi- 
eant species and plant groups. To avoid confu- 
sion between similar species, fresh samples 
were collected for groups of animals immedi- 
ately after the groups had vaeated an area. 
Bighorn sheep (Ovis canadensis) use steeper 
terrain on the northern range, and moose (Alces 
alces) are found at higher elevations; neither 
species was observed near the study sites. 


Shrub Utilization Rates 


Winter ungulate herbivory rates were sam- 
pled on the browsed transeets in late winter— 
spring before leaf emergence (usually late April) 
in 1963-1969, 1987, 1989, and 1990. Percent 
twig utilization was obtained from counts of al] 
browsed and unbrowsed twigs on each shrub 
located in the transeet. Diameters at basal point 
and browsing point were measured on 20 ran- 
dom shoots on every fifth browsed shrub of 
ach species, and bite sizes were estimated fol- 
lowing Pitt and Selwab (1990). 


Trends in Big Sagebrush, 
1958-1990 


Five exclosures were erected in 1957 and 
three more in 1962 (n = 8 total). The exclosures 
were placed in sites representative of mixed 
big sagebrush/bunchgrass communities. Paired 
belt transects (each 1.5 m X 30.5 m = 46.5 m2) 
were permanently located inside and outside 
eight of the exclosures (one per exclosure treat- 
ment) in big sagebrush communities (Canfield 
1941, Parker 195+). Each matched pair of 
transects was as nearly comparable as possible 
in terms of slope, aspect, elevation, shrub 
speeies, and shrub cover (Barmore 1980, 
Houston 1982); nevertheless, differences might 
have occurred. Sampling of transects occurred 
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at the date of exclosure, whieh should reveal 
any initial site differenees. The transect for 
exclosure was selected arbitrarily. Heights and 
species of all shrubs found on the belt tran- 
sects were recorded in 1958, 1962, 1967, 
1974, 19S1, 1986, and 1990. Numbers of indi- 
vidual shrubs and any shrub seedlings were 
tallied. Aerial cover of all shrubs was mapped 
on graph paper, and shrub cover was later esti- 
mated using a grid (Barmore 1980, Houston 
1982). 


Detailed Site Comparisons 
in 1986 and 1987 


Shrub belt transects were not replicated at 
a site (n = 1 transect per treatment per loca- 
tion, 13 transects total) and were useful pri- 
marily for long-term trends and assessment of 
pretreatment conditions (Parker 1954). In 1986 
and 1987 more intensive and better replicated 
measurements (1 = 15 plots per treatment) 
were gathered; 15 cireular plots, each 1.7 m in 
radius (9.3 m7), were randomly located in big 
sagebrush stands both inside and outside six 
exelosures. The tallest height, widest diameter, 
and perpendicular diameter were recorded for 
each shrub within eaeh plot. The number of 
totally dead shrubs was recorded. The percent- 
age of dead material on partially live shrubs was 
estimated. At every fifth shrub of each species, 
lengths and diameters of 10 randomly sampled 
twigs were measured, and every vegetative twig 
and reproductive stalk were counted. A mini- 
mum of 100 twigs of each species from each site 
was colleeted, dried, and weighed. Canopy 
area for each individual shrub was estimated 
following Peek (1970) from the formula for the 
area of an ellipse: 


area = Cr) did?, 


where d! = largest diameter and d? = its per- 
pendicular diameter. All plots and long-term 
transects were located more than 25 m from 
exclosure fences to avoid the effects of snow- 
drifts or ungulate trails along the fences. 

Mean shrub height, largest erown area, shoot 
numbers and lengths, and total shrub cover 
were compared using a two-way ANOVA, with 
browsing and exclosure loeation as treatments. 
The six exelosures should be considered repli- 
cations of one treatment (browsing) with 5 d.f. 
used to test for differences among the treat- 
ment and replieations. Replieations at a location 
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included the 15 random plots in each treat- 
ment (15 X 2 X & locations = 240 plots of 9.3 
m?). Nonparametric procedures were used for 
percent twig utilization comparisons between 
the 1960s and 1980s and for other data that 
were nonnormal or with unequal variances (F- 
max tests; Sokal and Rohlf 1981). Frequency 
distributions of shrubs in browsed and un- 
browsed plots were compared using the 
Kolomogorov-Smirmov test (Zar 1974). All dif- 
ferences discussed are significant at the P < 
.O5 level unless otherwise indicated. 

Aboveground biomass production of shrubs 
was estimated from the numbers of shrubs per 
plot times the average number of reproductive 
and vegetative shoots per plant times the aver- 
age dry weight of shoots. Regression equations 
for dry weight of shoots (independent vari- 
able) were calculated following MacCracken 
and Viereck (1990) from diameters at base (de- 
pendent variable) and length (dependent vari- 
able) of shoots. Separate regression equations 
were calculated for reproductive and vegetative 
shoots of big sagebrush, green rabbitbrush, 
and rubber rabbitbrush. The regressions on dry 
weight were applied to the sample of all twig 
diameters and lengths to estimate average 
twig biomass. 


RESULTS 
Ungulate Densities and Diets 


Ungulate densities approximately doubled 
during the study period on the low-elevation 
stratum (Table 1), whereas they approximately 
tripled on the high-elevation stratum during 
the same period (Table 1). Pronghorn consumed 
81% shrubs in their diet, followed by mule 
deer 50%, elk 8%, and bison 1% (Table 2). 
Pronghorn diets were 49% big sagebrush, 
mule deer consumed 23%, and elk diets were 
only 4%. The higher combined ungulate den- 
sities and the presence of pronghorn and mule 
deer, both of which eat more big sagebrush, 
suggest that ungulate herbivory on big sage- 
brush will be greater on the low-elevation 
stratum study sites. Since pronghorn consume 
12x more big sagebrush in their diets than elk 
and 2x more than mule deer, and since prong- 
horns also spend summers in the low-elevation 
stratum, we suspect pronghorn were the most 
important herbivore on big sagebrush on the 
low-elevation stratum. 


UNGULATE BROWSING IN YELLOWSTONE PARK 205 


TABLE 1. Average minimum densities of elk, pronghorn, 
mule deer, and bison near the lower- and higher-elevation 
exclosures on the northern winter range of Yellowstone 
National Park. Reported densities are based upon actual 
counts from fixed-wing aircraft (Houston 1982, Meagher 
1989, Singer 1991) and are uncorrected for visibility bias. 


Ungulate density (no./km2) 


Lower Higher 

Ungulate exclosures exclosures 
1965-19658 

Elk 6 6 

Pronghorn 3 0 

Mule deer 4 0 

Bison 0 I 
Total tle i 


1985-1985 


Elk 16-19 16-19 
Pronghorn Fale) 0 
Mule deer 2 0 
Bison tr 2 
Total 25-31 18-21 


Shrub Utilization Rates 


Big sagebrush utilization rates were consis- 
tently high (87%) and did not differ between 
1963-1969 and 1985-1988 at the low-eleva- 
tion stratum study sites dominated by the more 
palatable (to pronghorn) Wyoming big sage- 
brush (Table 3, Mann-Whitney U tests, P > 
.05). Pronghorn and elk reductions during 
1962-1967 apparently did not result in any 
decrease in percent leader use of Wyoming 
big sagebrush on the low stratum. Green rabbit- 
brush was also used heavily at the low-stratum 
sites where deer and pronghorn occurred (Table 
3). Utilization rates of big sagebrush at the 
higher sites dominated by mostly mountain big 
sagebrush, however, increased about sixfold 
after ungulates increased threefold (Table 3, P 
< .05). Use of green rabbitbrush did not in- 
crease significantly at the high stratum during 
this period of ungulate increase. Percent leader 
use of big sagebrush at the lower-elevation 
sites averaged 87%, but leader use averaged 
only 11% at the higher sites. Bite sizes averaged 
73% of vegetative shoots and 83% of repro- 
ductive shoots (n = 180 measured diameters 
of browsed shoots and 540 unbrowsed vegeta- 
tive and reproductive shoots). Consumption of 
annual aboveground biomass of big sagebrush 
by ungulates averaged about 68% at the low- 
elevation stratum sites and 9% at higher-ele- 
vation stratum sites. 
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TABLE 2. Mean percent of shrubs in winter die ts of four ace oy s on Yellowstone's northern winter range, 1985-1958, 
determine by microhistological analysis of feces (¥ + SE) 
( fou (no. 

aggregate Big Rabhit- Eurotia Fringed Total 

samples) sagebrush brush# lanata sage shrubsb 

sg SE av sy) x SE v BE T Sie 

Ly = Se ee eee Se ee ee ee ee 
Elk (28) 3.5 orl 1.6 2 0.4 1.0 1.0 20 10 4.1 
Bison (25) 0.1 0.] 0.0 0.0 0.0 0.0 0.0 0.0 1.3 1.6 
Mule deer (21) ie ee Vo.) 52 4.6 0.7 1.0 17.9 16.7 49.7 20.9 
Pronghorn (20) 48.7 15.0 ae Oe 3.5 6.0 18.5 oe $0.5 boat 





aRabbitbrush inchides Chrysothamnus nauseosus and C. visidiflorus. 
b}otal shrubs also includes Populus spp.. Salix spp., and Atriplex spp. 


Big Sagebrush Trends in Densities, 
Heights, and Cover, 1958-1990 


Big sagebrush on belt transects in the lower 
stratum differed at the time of exclosure in 
1958-1962. Densities were similar, but average 
heights were 50% and cover about 60% those 
values on transects selected for exclosure (Table 
4). Big sagebrush densities, heights, and cover, 
however, were similar between browsed and 
unbrowsed transects at the time of exclosure 
on the high-elevation stratum. 

Apparently, ungulates were suppressing 
Wyoming big sagebrush on the low-elevation 
stratum during the study period. Wyoming big 
sagebrush densities decreased 43% and big 
sagebrush cover decreased 29% on the low- 
stratum browsed site over the 3l-year period 
(Table 4). Density and cover of Wyoming big 
sagebrush increased dramatically (850% and 
§30%, respectively) in the unbrowsed sites of 
the low stratum. Big sagebrush individuals were 
taller on unbrowsed sites (Table 4). 

Herbivory effects were less on the high- 
elevation stratum study sites, and all trends 
were similar for browsed and unbrowsed sites. 
Mountain and basin big sagebrush density 
declined, and canopy cover increased on both 
browsed and unbrowsed belt transeets, 
1958-1990 (Table 4). Wyoming and basin big 
sagebrush density declined 39%, but cover 
increased 39% on browsed sites over the 31 
years. Heights of big sagebrush inereased on 
both browsed and unbrowsed sites, but more 
on unbrowsed sites (Table 4). 


Detailed Site Comparisons of Densities, 
Cover, and Biomass Production in 
1986 and 1987 
Densities of big sagebrush (F = 50.9), total 
canopy cover of big sagebrush (F = 8.1), indi- 
vidual shrub crown area (F = 22.5), and heights 


of big sagebrush (F = 79.8, P < .05) differed 
between a much larger sample of browsed (n 
= 15) and unbrowsed plots (n = 15 per loea- 
tion, n = 180 total) sampled in 1986 and 1987. 
In each case, however, location was also signif- 
icant, and the interaction between location and 
browsing was significant. For example, sage- 
brush individuals were 59% taller on unbrowsed 
plots at six exclosure sites, but at the Blacktail 
exclosures sagebrush plants were taller on 
browsed plots. Heights of big sagebrush, 
green rabbitbrush, and horsebrush increased 
with elevation in both treatments. As a conse- 
quence of this exploratory analysis and signifi- 
cant interactions with Joeation, our division of 
plots into a high and low strata appeared justi- 
fied, and we analyzed data from the lower and 
higher study sites separately in all subsequent 
analyses. 

Big sagebrush individuals were shorter and 
erowns smaller in browsed versus unbrowsed 


TABLE 3. Percent of twigs browsed in big sagebrush 
communities on Yellowstone’s northern range. Total un- 
gulate numbers increased twofold at the lower exclosures 
and threefold at the higher exclosures between 1963-1969 
and 1985-1988. The same transects of 46.5 m? each (2 = 5) 
were sampled both periods; only these five browsed tran- 
sects were sampled 1963-1969. 








Location 1963-1969 1985-1988 
Shrub species 
(n = transects) we SE a Si 
Low elevation (n = 2) 
Big sagebrush 8.0 4.2 86.8 T2 
Green rabbitbrush@ m1 * 105 
Spiny hopsage® 14.8 4.3 
High elevations (n = 3) 
Big sagebrush 1.9 0.8 11.6 Se 
Green rabbitbrush Gea Be) 8.9 2.9 
Horsebrush 46.6 IL-4 





“Only big sagebrush utilization was sampled 1963-1969, and green rabbitbrush 
at only the higher exclosnres. 
*P < .05, according to Mann-Whitney U tests. 
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TABLE 4. Changes in density, heights, and canopy cover of individual big sagebrush shrubs between time of exclosure 
placement in 1958 and 1990 on permanently marked 46.5 m2 shrub transects, Yellowstone's northern winter range. 


Heights Canopy cover 
Density of shrubs (cm) (m2/46.5 m2) 
1958-1962 1990 1958-1962 1990 1958-1962 1990 

Treatment ae SE x Si x SE <i SE x SE < SE 
Lower exclosures# 

Browsed 2) HH || 2 0.7 0.3 0.5 0.1 19 f 16 4 

Unbrowsed 25 15 103 27 1.8 0.6 16.7 0.9 28 wi 50 g 
Higher exclosures 

Browsed 67 1S 4] ee 1.9 0.9 5.9 1.6 12 2 42 7 

Unbrowsed GL 34 Sige: 6.2 1.6 0.2 S.6 0.9 10 2 62 9 


4Big sagebrush subspecies in these transects, mostly Af. teyomingensis, are apparently highly palatable to pronghorns. 


bBig sagebrush subspecies include mostly A. t. vaseyana. 


exclosure sites at low elevations (F = 29.8, 
14.3, respectively), but there was no difference 
in heights or crown sizes due to browsing at 
the high-elevation sites (P > .05, Table 5). 
Horsebrush was shorter and crowns were 
smaller on browsed and unbrowsed exclosure 
sites at the higher elevations only (F = 14.5, 
4.6, Table 5). Common rabbitbrush was short- 
er on browsed plots at the lower elevations, 
but it was taller on browsed plots at the higher- 
elevation exclosure sites (Table 5). Density of 
Wyoming, big sagebrush was less on browsed 
versus unbrowsed plots at the lower exclo- 
sures (F = 14.7), but there was no effect of 
browsing at the higher exclosures (Table 6). 
No difference in the number of dead big sage- 
brush individuals was observed between 
browsed and unbrowsed plots at either eleva- 
tion category (P > .05). Twenty-two times 
more seedlings of the year were observed on 
browsed than unbrowsed plots at higher ele- 
vations (F = 2.7, Table 6). 

Big sagebrush contributed 82-99% of annual 
aboveground shrub production in these shrub 
communities. Browsing did not consistently 
influence the production of big sagebrush or 
green rabbitbrush at higher exclosure sites, 
but browsed rubber rabbitbrush produced 
less biomass at higher-elevation sites (Table 6). 
Both Wyoming big sagebrush and rubber rabbit- 
brush produced much less aboveground bio- 
mass on browsed sites on the low study sites 
(Table 6). 

There was no influence from browsing on 
the number of vegetative or reproductive shoots 
per shrub for big sagebrush or green rabbit- 
brush. Reproductive shoots averaged 42% 
longer (Friedman test, y,? = 38, n = 6 locations, 
P < .05), and vegetative shoots averaged 45% 


longer on browsed versus unbrowsed big 
sagebrush (Friedman test, ¥,2 = 42, n = 6 loca- 
tions, P < .05, Table 7). There was no effect of 
browsing on length of reproductive shoots of 
green rabbitbrush (P > .05). 


DISCUSSION 


Other studies indicate mountain big sage- 
brush is preferred and eaten at a higher rate 
by mule deer and elk than Wyoming big sage- 
brush, while basin big sagebrush is the least 
preferred (Sheehy and Winward 1981, Welch 
et al. 1981, Personius et al. 1987). Our observa- 
tions initially appear in constrast with this gener- 
alization; we observed 70% more winter utili- 
zation on Wyoming big sagebrush than moun- 
tain big sagebrush. Too few basin big sage- 
brush occurred on the study sites to draw any 
conclusions. Our data do not constitute a pala- 
tability test, however, in that mountain and 
Wyoming subspecies did not occur at the same 
study sites. We suspect pronghorn were the 
primary herbivore on Wyoming big sagebrush 
in lower study sites; pronghorn find the Wyo- 
ming subspecies highly palatable (Beetle 1960, 
Beetle and Johnson 1982), and that subspecies 
was more available to all ungulates due to 
shallow snows and more winds in the low stra- 
tum. Ungulate preference for big sagebrush 
subspecies also varies between locales (Welch 
et al. 1981, McArthur and Welch 1982); for 
example, Dietz and Nagy (1976) found 
Wyoming big sagebrush was preferred by mule 
deer in Colorado. 

Mountain and basin big sagebrush seedling 
germination, establishment, and survival were 
apparently enhanced by browsing and ungu- 
late grazing (possibly due to secondary effects 
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TvBLe 5. Individual shrub crown and heights of shrubs in browsed and unbrowsed (protected) sites on Yellowstone’s 


northern elk winter range. Samples were drawn from (a = 


unbrowsed sites in 1986 and 1987. 


: ae 
Crown area (cin) 


180) plots of 9.3 m2? cach located randomly in browsed and 


Heights (em) 





Shrub species Unbrowsed Browsed Unbrowsed Browsed 
Location x 1) Be SE ay SE BY ro) 
Big sagebrush# 
Lower (Artrwy) 675 §5 347 66* 50 2 oF +* 
Higher 79S 90 2-4 ay) 79 3 fall z 
Horsebrush 
Lower 45 7 Cif 13 [iD ] ipl Z 
ERietier BE 222 ra 93* 63 10 27 3% 
Common rabbitbrush 
Lower 254 109 278 $5 7S § 43 on 
[igher 196 42 SSI 990* oY 3) 50 1] 
Green rabbitbrush 
Lower 76 13 104 oD 36 1] 28 3 
Higher 742 70 392 42 TY a D3 pea 


*Siznificant difference between grazed and control means using ANOVA, P < .05. 


aBig sagebrush subspecies included lower exclosures—A. t. teyomingensis only, higher exclosures—mixed populations of A. t. tridentata and A. t. vaseyana, but 


nearly all A. t. vaseyana. 


TABLE 6. Estimated annual production (g/m?) of the most common shrubs in browsed and unbrowsed big sagebrush 
communities at six exclosures on Yellowstone's northern winter range (n = 15 plots each in both browsed and 
unbrowsed treatments at each site). Wyoming big sagebrush is found only at the lower-stratum exclosures, and mixed 
populations of nearly all mountain with some basin big sagebrush at the higher exclosures. 


Lower e] evations 


Higher elevations 


Unbrowsed Browsed Unbrowsed Browsed 
Exclosure loeation a SE By SE a SE x SE 
Estimated biomass 
(g/m?) 
Big sagebrush 1S 1.9 73.6 72.6 
Green rabbitbrush Onl 125 ie) 5.8 
Common rabbitbrush 0.6 0.3 5.9 a 
No. big sagebrush 
individuals/9.3 m2 
No. alive 16 2 2 1+ 13 2 Vs: ) 
No. dead 13 Ot 0.6 0.5 a6 0.5 5.9 1.9 
No. seedlings 0.8 0.3 02 0.1 Ow 0.1 4.4 es 


*P < .05 according to t tests. No tests were conducted on biomass since it was estimated from a product of no. of plants X average no. of shoots < average 
weight of shoots. Tests were conducted on each of those parameters separately, however (see text and Table 7), suggesting statistically significant differences at 


the lower elevations. 


such as reductions of herbaceous vegetation 
competition and ungulate hoof action) at the 
higher winter range, but the opposite trend 
was observed on lower sites. MeArthur et al. 
(1988) also observed more big sagebrush seed- 
lings on a site browsed by mule deer in winter 
than on an unbrowsed site. The physical act of 
ungulate grazing, with its accompanying hoof 
action, greater soil disturbanee, more bare 
ground, and less standing dead vegetation and 
litter, may provide conditions more suitable to 


big sagebrush germination. Big sagebrush indi- 
viduals are smaller on browsed sites, which 
may also benefit establishment and survival of 
seedlings due to reduced competition for 
light, soil moisture, and other resources. 
Ungulate herbivory suppressed big sage- 
brush on the lower-elevation sites, where almost 
no recruitment of Wyoming big sagebrush 
occurred on browsed sites; apparently few sced- 
lings survive the intense browsing. Wyoming 
big sagebrush reproduces more suecessfully 
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TABLE 7. Numbers and lengths of reproductive and vegetative stalks on shrubs in browsed and unbrowsed plots in big 
sagebrush communities on Yellowstone's northern range (n = no. shrubs), 





No. reproductive flowers/shrub 


Length (cm) of flower stalks 


Unbrowsed Browsed Unbrowsed Browsed 

Species 

Location x SE x SE x SE z yo 
Big sagebrush 

Lower (Wyoming spp. only) 15 6 13 10 4 0.2 ts) 0.5%" 

Higher (basin and 

mountain spp.) 2) 7 17 6 14 0.7 18 0.6** 

Green rabbitbrush 

Higher 10 3 16 10 U1 O04 12 0.4 


No. vegetative shoots/shrub 


Big sagebrush 





Lower 99 a3 8S 18 2 0.1 a 0.4 

Higher 83 3) a9 12 5 0.5 7 0.5** 
Green rabbitbrush 

Higher 2g ji 52 20 6 0.2 8 0.4** 
JUS (Sy. 


**P < 01. Differences between numbers in browsed and unbrowsed plots were tested with t tests and Jengths with Mann-Whitney U tests. 


than the other subspecies on xeric sites (Welch 
and Jacobson 1988), and the xeric, sodic clay 
soils of the low stratum are clearly more suit- 
able to Wyoming big sagebrush. At the high lev- 
els of ungulate herbivory we observed (rough- 
ly 68% biomass removal), the Wyoming sub- 
species is presently suppressed by ungulates. 

The ability of Wyoming big sagebrush to 
recover from herbivory is less than for moun- 
tain and basin big sagebrush. Wyoming big 
sagebrush is shorter (individuals often do not 
exceed 0.3 m), seedling growth rates are lower, 
and current annual growth is less than for the 
other two subspecies (McArthur and Welch 
1982, Booth et al. 1990). The approximately 
66% decline in numbers of mule deer using the 
lower stratum within the park over the past 
two decades may be due to the localized Wyo- 
ming big sagebrush decline. Pronghorn did not 
decline in the lower stratum during the same 
period, but pronghorn, unlike mule deer, were 
artificially reduced well below carrying capac- 
ity levels during the 1960s (Houston 1982), 
and they may still be recovering from the 
reductions. 

Increases in height and cover of big sage- 
brush are reported after protection from ungu- 
lates. Robertson et al. (1970) reported big sage- 
brush cover increased 76% after 30 vears of 
protection from browsing, although mean 
heights declined 12%. Heights and crown sizes 
were similar, but live cover by big sagebrush 
was greater on unbrowsed sites on a mule deer 


winter range, primarily due to a greater die- 
back of browsed big sagebrush (McArthur et al. 
1988). Average crown dieback was 64% in the 
browsed area and 17% in the unbrowsed area 
(McArthur et al. 1988). Mule deer use was heavy 
(370 deer-use days/ha), and dieback of big 
sagebrush occurred after two successive win- 
ters of heavy snowfall (McArthur et al. 1988). 
Browsing by native ungulates stimulated 
seedstalks and leaves of big sagebrush and 
leaves of green rabbitbrush on the study sites. 
Stagnation of shrubs occurred inside big game 
exclosures after only two years of exclosure— 
nonuse of big sagebrush resulted in an aver- 
age 36% reduction in biomass production over 
clipped plants (Tueller and Tower 1979). 
Numbers of sprouts of green rabbitbrush were 
similarly increased by clipping (30% herbage 
removal), and new growth was longer, leaves were 
larger, and leaves remained green for one 
month longer (Willard and McKell 1978). On 
the other hand, browsing of more than 80% of 
the leaders of mountain big sagebrush by 
mule deer resulted in a reduction of 50-93% 
in total number of seedstalks per plant and a 
reduction of 0-53% in length of seedstalks 
(Wagstaff and Welch 1991). Grazed grasses 
on the northern Yellowstone winter range 
have higher protein levels (Coughenour 1991), 
and grazing stimulates aboveground growth 
of grasses (Frank and McNaughton 1993). 
Increased vigor in new growth of browsed 
shrubs on the Yellowstone northern winter 
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range is consistent with these observations of 
grasses, and shrub vigor may be the result of 
increased rates of nutrient cycling due to 
ungulate defecation and urination (McNaughton 
1979). In addition, plant competition is reduced 
and water availability may be increased on 
browsed sites on the northern winter range 
due to smaller crown sizes and fewer transpir- 
ing tissues for individual shrubs. 

Historic mean duration between fires was 
25 years on the Yellowstone northern winter 
range, but due to fewer fire starts and active 
fire suppression, no significant burning of 
grasslands occurred between 1870 and 1988 
(Houston 1973, Romme and Despain 1989). Big 
sagebrush communities had not yet achieved 
climax postfire state on the northern winter 
range as indicated by increases in heights and 
cover of both browsed and unbrowsed big 
sagebrush individuals between 1958 and 1990. 
Browsed big sagebrush communities on higher- 
elevation ranges were replacing themselves; 
many successful seedlings and small individuals 
were observed on browsed versus unbrowsed 
sites. Lomasson (1948) observed almost no 
reproduction for 40 years in a stand of big 
sagebrush, but then reproduction increased as 
the original population began dying. Average 
life span of big sagebrush is 53-72 vears, and 
in amature, undisturbed stand, most big sage- 
brush individuals are in the 55-59-year age 
class (Roughton 1972). Sagebrush recovery 
following fire varies from a few years to 30 
years depending upon environmental condi- 
tions for reestablishment (Sneva 1972, Harniss 
and Murray 1973). [f most big sagebrush com- 
munities we studied on the northern winter 
range last burned in the 1840-1890 period 
(Houston 1973), then most big sagebrush pop- 
ulations should have approached senescence 
and population turnover at the time of the 
1986-87 sampling. 

Ungulate herbivory levels on the lower 
study stratum restricted growth, establish- 
ment, and survival of big sagebrush at the 
time of this investigation, although browsed 
big sagebrush communities were stable or in- 
creasing at the higher elevation. Suppression 
of growth and reproduction of plants by in- 
creasing native ungulates can result in a new, 
altered plant-ungulate equilibrium (Sinclair 
1977, Caughley 1981). If unnatural (human- 
caused) concentrations of ungulates cause 
plant alterations, the situation is not accept- 
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able under NPS policy (U.S. Department of the 
Interior 1988). Houston (1982) concluded 
densities of ungulates in the BLA were unnat- 
ural and artificially high due to animal avoid- 
ance of hunting outside the park. If so, some 
form of ungulate management—control, en- 
couragement of migrations—is justified on the 
BLA. Elk and pronghorn reductions in the 
1960s, however, did not reduce percent leader 
use or improve the declining status of big 
sagebrush in the BLA. Either effective ungu- 
late densities remained the same near the big 
sagebrush study sites, the ungulate reductions 
did not go on long enough, or high preference 
for the Wyoming subspecies by pronghorn 
maintained high levels of herbivorv in the area 
during the control period. 

We caution that we were unable to calculate 
appropriate or recommended ungulate herbi- 
vory levels or ungulate densities for the north- 
em winter range. Our data included two dichot- 
omous periods in ungulate management. The 
first period of our study, 1958-1968, was clear- 
ly a period of ungulate underpopulation dur- 
ing which time elk, bison, and pronghorn 
were controlled far below ecological carrying 
capacity (ECC) densities (Barmore 1980, 
Houston 1982, Boyce 1993, Mack and Singer 
1993, Singer and Norland 1995). The second 
period of our investigations, 1986-1988, likely 
was a time of ungulate densities in excess of 
natural conditions, at least for elk and bison. 
This statement is not based upon any compar- 
isons to control conditions (no similar ecosystem 
exists with wolves and nondisrupted migra- 
tions for a comparison), but upon the conclu- 
sions of Houston (1982) that elk concentra- 
tions were unnaturally high in the low-eleva- 
tion BLA stratum, and computer predictions 
that elk and bison would number 8-25% less 
following wolf restoration (Garton et al. 1990, 
Boyce 1993, Mack and Singer 1993). Pronghorn 
densities in relation to ECC are unknown— 
one author feels coyotes (Canis latrans) are 
suppressing pronghorn on the northern Yellow- 
stone winter range and that, following wolf 
restoration, coyotes will decline and prong- 
horn will further increase (Berger 1991). Wolf 
restoration occurred on the study area during 
the winter of 1994-95, providing an opportunity 
to test the effects of wolves upon ungulate- 
plant interactions in the Yellowstone ecosys- 


tem (Cook 1993). 
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